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Abstract 

We discuss the motion of extended objects in a spacetime by considering a gravitational field 
created by these objects. We define multipole moments of the objects as a classification by Lie 
group SO(3). Then, we construct an energy-momentum tensor for the objects and derive equations 
of motion from it. As a result, we reproduce the Papapetrou equations for a spinning particle. 
Furthermore, we will show that we can obtain more simple equations than the Papapetrou equations 
by changing the center-of-mass. 
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I. INTRODUCTION 



The motion of celestial bodies has attracted people's interests from the ancient era. The 
ancient people used to obtain the important information on the seasons or the time from it. 
Even now, the motion of celestial bodies is one of the important foundations of the modern 
science. The modern physics begins with the Kepler's law, the Newton's law and the law 
of universal gravitation. Our understanding about the celestial motion was extended by the 
Newton's theory. The interest in the celestial motion did not change even when Einstein 
proposed the general relativity. Indeed, one of the classical tests of the general relativity is 
the account for the anomalous perihelion shift of Mercury [lj]. At present, the knowledge 
of the celestial motion is increasing its importance as the devices to probe our universe. 
For instance, the general relativity predicts the existence of gravitational waves, and many 
researchers in the world endeavor to catch them. One of the main targets in the detection 
of gravitational waves is a binary system which consists of celestial bodies like neutron stars 
or black holes. Those researchers attempt to detect the gravitational waves generated by 
the orbital motion of the binary. Therefore, it is still important to investigate the motion of 
extended bodies like celestial objects in the gravitational field DQD 

Now, it is recognized that there are two flows in the studies on movement of an extended 
object in the general theory of relativity. One of them starts from the study of Einstein, 
Infeld and Hoffmann ^|. They regarded the existence of object as singular points of the 
gravitational field, and discuss the motion of objects in a weak field approximation. The 
advantage of their method is that it is unnecessary to treat the objects as test particles 
because Einstein et.al. solved the field equations in an approximation. Their study has led 
to the later development of the post Newtonian theory. Another flow regards the motion 
of objects as motion of a test particle on a background spacetime. In this standpoint, Pa- 
papetrou derived equations of motion for a spinning particle moving in a curved spacetime 
by considering the Taylor expansion, that is, multipole expansion of the divergence of an 
energy-momentum tensor at the position of the particle The equations derived by him 
are known as the Papapetrou equations. The Papapetrou's method was not covariant in in- 
iegrating tensors, but Dixon developed a covariant method by using the bitensor formulation 
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In addition, there exists another interesting way to derive the Papapetrou equations. 
Tulczyjew derived the Papapetrou equations in a more simple and covariant way by using 
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the Dirac's delta function |8[. More recently, Anandan, Dadhich and Singh developed an 
approach with the action principle |9]. They consider the Taylor expansion of an action for 
a particle and derive equations of motion. In this article, our standpoint is the one of latter 
flow. 

Now, there are two approaches to treat the multipole properties of objects. One of them 
is based on the distribution of the matter, and another one is based on the gravitational field 
created by the objects. In consideration on the motion of objects, it is indeed intuitive to 
focus on the distribution of matter in deriving equations of motion as Papapetrou and Dixon 
did. However, when we discuss the motion in the gravitation field, it is more natural to be 
based on the field created by the object than to be based on the distribution of the matter. 
In addition, the treatment of multipoles in previous researches was not exactly based on 
50(3), though the concept of multipole moments should be defined by a classification by 
the Lie group 50(3). 

In this article, our purpose is to discuss the motion of objects on a curved spacetime 
based on the above viewpoint. Our argument will be presented as follows. In the section HU 
we classify solutions of the field equations and define multipole moments based on 50(3) in 
the context of the special relativity. Then, we recognize an object as a point singularity of 
the field and construct an energy-momentum tensor of the object. In section IIII1 then we 
turn to the general relativity and derive equations of motion from the energy-momentum 
tensor. Then, when we examine up to dipole moment, we show that our equations can be 
equivalent to the Papapetrou equations in adopting a certain definition of the center-of-mass. 
Moreover, we propose more simple equations of motion for the spinning particle. In section 
IIV1 we summarize our results. 

Notations used in this article are the following. The signature of metrics is (— , +, +, +). 
We denote metrics of a Minkowsky spacetime and a curved spacetime by r]^ u and g^ u , 
respectively. The Roman indices run over 1 to 3 and the Greek indices run over to 3. The 
definition of curvature is R! 1 VOL p = - + T^^V 7 ^ - ^T va and R^ u = R 7 ^. 

We use convenient symbols d a ^...^ = <9 a c) 7 ■ ■ ■ d 1 for simplicity of writing. We occasionally use 
the multi-index notation for capital Roman indices; in the expression A^B L , the upper-case 
letter "L" indicates an ordered set of lower-case letter "/" indices, i.e., "I" is a mere number 
of the indeces. For example, A^B L = Aj tl ...j tl B kl "' kl . The signature of perfect asymmetry 
symbol e's is defined by e 123 = 6123 = +1 for 3-dimensional and e 0123 = — 60123 = +1 for 
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4-dimensional. 



II. MULTIPOLE PARTICLE IN SPECIAL RELATIVITY 

In this section, we formulate an energy-momentum tensor which is corresponding to a 
multipole particle in the inertial motion on the Minkowsky spacetime. Here, we suppose 
that the particle can be a source of gravitation but its motion receives no influence from the 
gravitational field. The interaction between the particle and the gravitational field will be 
introduced in the next section. 



A. Isolated solution of linearized Einstein equation 

It is well known that the Einstein equation can be written in a Klein-Gordon type equation 
as following, 

UhT = (-<?) (TV + t^) + (hT^hT - h^^hr, p ) , (2.1) 

where □ is the d'Alembertian by the Lorenz metric rf v = diag(— 1, 1, 1, 1), h^ u is defined 
as = rf* — yj—ggv, and t^L * s the Landau-Lifshitz's pseudo-tensor. Moreover, the de 
Donder gauge condition h^ u jU = is imposed on W as a coordinate condition. If we ignore 
the nonlinear terms in the equation (|2.1j) . we obtain the linearized Einstein equation 

1 fiirC 

UhT = — —r-TV. (2.2) 

The concept of multipole moment for isolated solutions of the equation (j2.2j) is given by 
a classification based on the representation of Lie group SO (3) Q- When the elements of 
SO (3) act on a spatial hypersurface which is spanned by the vectors d x , d y and d z , the h 00 
component behaves as a scalar, the h° l as a vector and the as a 2-rank symmetric tensor, 
respectively. 

Now, let us investigate the solution of the equation (|2.2j) under the assumption that the 
matter exists locally in the neighborhood of the origin. Namely, we assume that the matter 
is confined in the inside of a (4-dimensional) tube with the radius R and there is no matter 
in the outside of the tube (see the figure [TJ. Additionally, we adopt the boundary condition 
that hy- v — > as r — ► oo. Under these assumptions, the exterior solutions of the equation 
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FIG. 1: The tube and the distribution of the matter 



(|2.2|) is given as follows, 
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where we introduce a retarded time t re t = i — r. The M L etc. in ()2.3)1 are the integral 
constants and they behave as 3-dimensional tensors. It should also be noted that all the 3- 
dimensional tensors M L etc. are perfectly symmetric and traceless: M kl "' kl = M^ kl '" k ^ and 
M % i L bij = etc. We call the M L 's the gravitational multipole moments. Each gravitational 
multipole moment is a element of an irreducible representation space of SO (3) because of 
its symmetry. 

On the other hand, the interior solution is determined by the distribution of the matter 
in the tube. By jointing the interior and the exterior solutions at the surface of the tube, we 
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can determine the gravitational multipole moments M L 's. The important point is that any 
object in the exterior gravitational field of the tube cannot distinguish any interior structure 
if there is no difference in the M L 's. In other words, the only M L, s prescribe all the external 
gravitational properties of the matter in the tube. According to this view, we can freely 
deform the distribution of the matter in the inside of the tube as long as the deformation 
does not change M L, s. Among such deformations, the simplest one is contraction of the 
matter into a point. This can be achieved by extrapolating the exterior solution into the 
inside of the tube (this is equivalent to R — > 0). By this deformation, the extended matter 
is replaced with a particle that is represented by a singularity of the gravitational field. We 
call this particle as a multipole particle. Thus, our objective now is to investigate the motion 
of this multipole particle. However, before proceeding further, we examine the properties of 
the solution ()2.3|) . It is because the solution ()2.3|) has some gauge freedom and we have not 
used any gauge condition in (|2.3|) until now. 

Firstly, we obtain the following relations among the tensors from the de Donder condition, 



M L + P L + Q L = 0, 



(2.4a) 



P L + B L + -D L 
2 



(2.4b) 



S L + -C L + —E L 
2 2 



(2.4c) 



Q L + A L + -D L + F L = 0, 



(2.4d) 



where dots over tensors show the derivatives by t. Next, since the general form of ()2.3|) still 
includes some gauge freedom, we can further restrict the form of the solution by another 
gauge transformation. Let us consider the following gauge transformation 
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After this transformation, we become aware that the tensors A L , D L , E L and F L can be set 
to 0. Consequently, we obtain the following expression for a general isolated solution of the 
linearized Einstein equation, of which all gauge freedom is fixed: 
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where we define new symbols J M ' L_1 = e iam s mL ~ 1 . By these equations, we define the multi- 
pole moments of the field created by the multipole particle. We call the M L mass multipole 
moments (or mass multipoles), and the J iJ ' L spin multipole moments (or spin multipoles). 
Note that the mass multipoles begin with mass monopole, but the spin multipoles begin with 
spin dipole. The mass dipole moment, especially, decides the relation between the center- 
of-field and the mean position of the objects. If there is any necessity to distinguish S L and 
J l *' L , we refer them as original spin multipoles and variant spin multipoles, respectively. 
However, we will use the variant spin multipoles in almost all the case. 



B. Multipole moment in oblique frame 

In the last subsection, we obtained the multipole expansion of the linearized gravitational 
field which is based on a representation of SO(3). However, in the definition of multipole 
moments, there is still some freedom of selection of 3-dimensional hyper-surface which 50(3) 
acts on. In other words, "on which 3-dimensional space we define the multipoles?" Accord- 
ingly, what we should consider here is how we can define multipole moments in an inclined 
3-dimensional space which is not orthogonal to the 4-velocity of the particle. 

In order to specify a spatial 3-dimensional subspace, we identify a timelike vector, namely 
v, which is orthogonal to this subspace. Simultaneously, we consider {e{\ (i = 1, 2, 3) as unit 
vectors which span the hyper-surface. Then, {v, e{\ forms an orthogonal tetrad. We call it 
(v, e)-tetrad frame. On the other hand, we denote the velocity of the particle by u and set 
\ — di (i — 1, 2, 3). We call the orthogonal tetrad {u, Aj} as (u, A)-tetrad frame. These two 
tetrad frames, namely (u, A)-frame and (v, e)-frame, are related by a Lorenz transformation 
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FIG. 2: Definition of frames. 
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where x M is the coordinate of (it, A) -frame and x^ is of (v, e) -frame. Furthermore, let a 
tetrad constructed by {u, e^} be referred as (u, e)-tetrad frame. This (it, e)-frame is the 
oblique frame that we should consider (Fig. EJ). 

We cannot apply the similar discussion in the last subsection directly to the oblique case 
because the non-orthogonal property makes the gauge transformation too complicated to 
find an appropriate gauge function However, we can give general isolated solutions in 
(u, e)-frame as below, which is almost the same expression ()2.3|) in (u, A)-frame: 
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Here h^ u is the components of h in (u, e)-frame and d a is the derivative along the vector e a . 
Furthermore, we introduce t Tet and f as the retarded time and "oblique" radius, respectively. 
If the Lorenz transformation between (u, A) and (v, e) is the boost in x-axis direction, namely 
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t = at + [3x, x = pi + ax, (a 2 — f3 2 = 1), then they become f = a/ (ax) 2 + y 2 + z 2 and 
t re t = t+/3x—r, where x, etc. is the coordinates in (u, e)-frame. However, concrete expression 
of t and r is not important. The existence of an one-to-one relation between t r et ; ? and £ rct , 
r is essential. 

We focus on the fact that both ()2.3|) and ()2.9|) are the general isolated solutions. We can 
define a one-to-one mapping between these solutions as following, 
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Given that this mapping is one-to-one, we can construct a representation of SO (3) on (it, e)- 
frame by using the representation p on (w, A)-frame as following, 



p = T- l pT. 



(2.12) 



According to this construction of the representation, we can consider that W" is also a 
gauge-fixed solution if T~ x [h) iv {t xeX S) — h^ u (t rct ) is gauge-fixed. Therefore, we obtain the 
general isolated and gauge-fixed solutions in (u, e)-frame under the representation ()2.12j) . as 
following, 
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We define multipoles in the oblique frame by these equations in the same way of the normal 
frame case ()2.7jl . The important point is that the multipoles in ()2.13|) are perpendicular to 
v and not u, though those in ()2.7|) are perpendicular to u. 



C. Energy-momentum tensor for multipole particle 



Here, we will derive an energy-momentum tensor for a multipole particle from the solution 
of linearized Einstein equation. The relation, 
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Here, □ is the d'Alembertian in (u, e)-frame and we use commutability between □ and d^. 
This energy-momentum tensor is in the comoving frame of the particle. By these equations 
()2.15|) . we define multipole moments of a particle. In other words, the multipoles of the 
particle is defined by the multipoles of the field created by the particle. In this definition, 
there is no direct relation between the multipoles and the distributions of the matter which 
constructs the object. We regard the point r = as the position of the particle, but this 
position is unrelated to the mean position of the object. It is the mass dipole moment to 
decide the relation between the loci of the particle and the object. We consider that the 
particle represents the object by fixing the mass dipole moment. Accordingly, we investigate 
the motion of the particle hereafter. In order to get the energy-momentum tensor in a general 
frame, we only have to transform the coordinates as 



(2.16) 
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After the transformation, we obtain the following expression, 
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where we introduce a covariant delta function V(x) = ^ = and use a relation, 

5 3 (f) = J V(x" - z»)dr. (2.18) 
The multipoles in the general frame are defined as follows, 

M pi-Pn = gPi . . . e ^M ir "S (2.19a) 

J^-Pl-Pn = e P e » . e Pl . . . ^njab^-i^ ( 219h j 

Noting that the 3-dimensional epsilon t %ik in (u, e)-frame becomes v\ e J^ in a general frame, 
and we obtain the following relation, 

jTO-Pn = Vx ^ =g a sS PPl ~- pn . (2.20) 

Because we have defined the multipoles on the inclined 3-dimensional subspace, the multi- 
poles must be orthogonal to the normal vector v of the subspace: 

v p M»- = 0, (2.21a) 

VpJ^o- = 0, (or VpS*- = 0), (2.21b) 

and JP p -P~- obeys 

v ^-p- = (2.22) 

from the asymmetry of e Xflua . The relation for the case of or is known as an spin 
supplemental condition. 
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III. MULTIPOLE PARTICLE IN GENERAL RELATIVITY 



In this section, we will discuss the energy-momentum tensor for the multipole particle 
and its motion on a curved spacetime. 



A. Introduce action between particle and gravitational field 

It is easy to obtain an energy-momentum tensor on a curved spacetime from that given 
in the last section. This can be achieved by the "colon-to-semicolon" rule. Interpreting 
the rule from the perspective of the gauge theory, it introduces the interaction between the 
particle and the gravitational field by the gauge principle. However, the way to apply the 
rule to our case is not unique because the multipoles etc. and the covariant derivatives do 
not commutate each other. This issue could be grasped as the identification of a self-field 
of the particle on the curved spacetime because the difference by the orders of them is a 
combination of the connections and their derivatives. Therefore, we consider that the order 
of the multipoles etc. in ()2.17|) identifies the self-field and we will not discuss this problem 
in detail in this article. Accordingly, we can give an energy-momentum tensor for a particle 
on curved spacetime as following, 

oo 



n=0 



OO p OO p 

- W rfrV Pl ... p „_ 1 (^M^' ^- 1 I?)-^ / drV Px ... Pn _ x (u v M^ ^v) 
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Here, the dots over the multipoles show the derivatives by the proper time r of the particle, 



— = vTV^ (3.2) 



D 
Dr 

and the multipoles must also satisfy the slicing conditions (|2.21|) and (J2.22J) . We can ob- 
viously obtain the well-known expression for a point particle from our energy-momentum 
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tensor (|3.1|) when we take the only mass monopole moment. 

The equation of motion of the particle can be derived from T^ u . v = that is the inte- 
grable condition of the Einstein equation. However, we give an algorithm using Tulczyjew's 
theorems instead of using it directly to obtain the equation of motion. 

Lemma 1 For an arbitrary point x in the spacetime and for an arbitrary tensor a?'" defined 
on the world line of a particle, the following equation holds, 

/+oo 
V a • • • V^Vp [a^-{r)uP{r)V{x - z{r))\ dr 
-oo 

dr, (3.3) 



V Q • • • Vp 

OO 



-±±V{x-z{t)) 

Vt 

where w M is J^-velocity of the particle whose position is denoted by z{r). 

Theorem 2 (Tulczyjew) Let n fll "' flkUl "' Um (s) be an arbitrary tensor on a world line and 
N vl '" Vm be a tensor which is constructed from n fll '" flk ' /1 '"' /rn as follows 

n r+oo 

N w~vm =J2 V w • • • V Mfc [n^~ m " Vm V{x - z{r))\ dr. (3.4) 
fc=o 

Then N Ul "' Um can be reduced to a canonical form 

n r+oo 

N u 1 -u m = J2 / ■ ■V, k [v fll - flkUl - Um (r)V(x - z(r))}dT, (3.5) 

k=o 

where v^-wi-v™ = v (ni-n k >vv m anc / u ^ v ^i-^i-^ m = q 

Remark that it is not necessary for uW'Wi" - ^ to be orthogonal to u M by the contraction with 
the indices z/'s of . Moreover, is not generally equal to n^ 1 '"^" 1 '""™ . 



Theorem 3 (Tulczyjew) If N Ul ' Um in theorem^ satisfies 

[ N^" Vm K vl ... Vm sTg d 4 x = (3.6) 
Jd 

for an arbitrary K Vl ... Vra and an arbitrary domain D, then all the coefficients i)^"^ v ^" Vm in 
the canonical form of N Vl ' Um must vanish. 

These theorems were given in Tulczyjew's paper. However, proofs of them were omitted or 
a little imperfect in his paper, we give them in the appendix of this article. Using these 
theorems, we propose an algorithm as follows: 

13 



1. Calculate the divergence of the Energy-momentum tensor. 

2. Using the lemma ^ and the commutation relation of the covariant derivatives, trans- 
form the result of step ^ into the form in the theorem EJ 

3. Using the theorem 01 derive the equation of motion from the coefficients of the canon- 
ical form. 

We can systematically derive the equations of motion by these steps. 



B. spinning particle 

Here, we consider a spinning particle as a non-trivial example. The spinning particle has 
angular momentum, and so we can obtain an energy-momentum tensor for it when we take 
n = 0, 1 in the equation (JXTJ): 

= J dr { (Mu p u v - M p u v - M u u^ V + V p (M p u p u v + u" J pp + u p J up ) V } . (3.7) 

In this stage, we still leave the mass dipole moment M p not to be fixed. What is to do is to 
process (13. 7|) according to the algorithm in the last subsection. For the simplicity, we set 

A"" = Mu"u v - M p u u - M v u p , (3.8a) 
B p V » _ m p u p u u + u u J pp + u^J up . (3.8b) 

Moreover, we decompose A up and B pvpj into the compositions parallel and perpendicular to 
u as follows, 

B p»h = b P^ _ u v b P m " _ u P b ^ + u p u "\f^ (3.9b) 
where we define the following variables, 

a vp = P\A xp , (3.10a) 

a p = u u A Ufl , (3.10b) 

b pvp = P p a P\B uXp , (3.10c) 

bf = P p a u x B aXp , (3.10d) 

b vp = u a P\B aXp , (3.10e) 

b p = u a u x B aXp , (3.10f) 
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and P^ v is the projection operator which is defined by P^ v = 5^ v + u^u v . After some 
calculation, we obtain the divergence of the energy-momentum tensor, 



T pv , v = J dr\ V U V P (b^V) + V, 



q ph + P\R- ( u a b p - b a / - b a /) ) V 
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Consequently, by using the theorem El the following equations are acquired: 



3H 



(3.12a) 



a PP + pp (u°b» - b\ p - b a /) = 0, 
Dr 



(3.12b) 



D 
~Dr~ 



-a" - u 



D 



bT - + [V„ Vj Qfo^ - u"b^j = 0. (3.12c) 

Let us substitute ()3.8|) and ()3.1U|) for ()3.12|) . Firstly, the equation ()3.12a|) becomes an identity 
= 0. Secondly, the equation ()3.12b|) turns into 



-M a u^ + — {-J ap + M a u^) 
Dt 



0. 



By multiplying u„ to this, we obtain 



J ap u p = 0. 



By multiplying P u tl to (J3.13|) and symmetrizing the indices v and p, we get 



P\P p a (M {a ii p) ) = 0, 



and by antisymmetrizing, we obtain 



D 



P\P P a {-J*" + M [rT U P] ) = 0. 

Dt 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



Thirdly, the equation ()3.12c|) with ()3.14|) changes to 

[Mu p -M p - (u a M a ) u^+^R p ia ^J aP -R p a ^J aP v?-R p ia ^M a v? = 0. (3.17) 
Because the relation i? M aj a 7 + R p p ia + R p -yaf3 = concludes 



D 
Dt 



R p aM J af3 v? = -\R p ^v?J a ^ (3.18) 
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we can reduce the equation f)3. 12c|) to 

(Mil* -M»- (u a M a ) u") + R^apu 1 (J al3 - M a u p ) = 0. (3.19) 

Summarizing the result, we obtain the equations (|3.14|) . (|3.15p . ()3.16|) and ()3.19|) as the 
equation of motion. 

Now, we left the mass-dipole M M not to be fixed until here. By choosing a suitable defini- 
tion of the center-of-mass, we can show that the equations are equivalent to the Papapetrou 
equations, 

(mil" + u a S a ^ + ^R" yaf3 u^S af3 = 0, (3.20a) 

gr> + S^UaU" ~ S^U^Il" = 0. (3.20b) 

To do this, we redefine the mass monopole and spin-dipole as follows, 

= 2J" V - M"u u + M u u" (3.21a) 



m = M+ [u a M a j (3.21b) 

It is obvious that (|3.1fi|) becomes to (|3.20b|) under the equation ()3.21aj) . On the other hand, 
the replacement ()3.21b|) makes the contents in the first parenthesis of (|3.19|) to 

MvP — M" — (u a M a ) u* 1 = mu" - P" a M a - (u a M a ) (3.22) 

Now, notice the fact that we can set 

P" c M a + (u a M u ) u" = S^Ua, (3.23) 

because of the freedom of M M . The mass-dipole M M is orthogonal to v, so that, M M has 
three degrees of freedom. On the other hand, since the equation ()3.23|) is orthogonal to u, 
()3.23|) means three independent equations. To sum up, the equation ()3.23|) is a set of three 
ordinary differential equations for the three independent components of M". According 
to the theorem on the existence of the solution of the ordinary differential equations, we 
can assure the existence of M M satisfying (|3.23j) . Consequently, we have showed that the 
equation ()3.19|) is equivalent to the Papapetrou equation (|3.20a|) . Moreover, we can also 
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show the energy-momentum tensor (|3.7|) is equal to the energy-momentum tensor given by 
Tulczyjew Q], 

r f / 1 DSw 1 DS up \ 1 1 

(3.24) 

under the replacement (j3.21ajl . ()3.21hjl and 

In the above discussion, we set M M to satisfy ()3.23j) . However, This is not only a unique 
choice. Rather, the equations become more simple by setting M M = as well as the case of 
Newton gravity. When we set M M in this way, we conclude the following equations 

M— - + R» aP u~<J a P = 0, (3.25a) 
Dt 



D 

= 0, (3.25b) 

Dt 
DM 



Dt 0, (3.25c) 

as the equation of motion of the spinning particle. The equation ()3.25a)l was known as 
the pole-dipole approximation of the equation of (|3.20a|) . but, it is no longer necessary to 
mention pole-dipole approximation. 



IV. SUMMARY 



In this article, we constructed an energy-momentum tensor for a multipole particle, by 
considering a gravitational field created by the object. The multipole moments of the object 
were defined by those of the gravitational field, and they were elements of the irreducible 
representation space of 5*0(3). In this definition of the multipoles, there were two kind 
of arbitrariness. One of them is a slicing condition that decides a spatial hyper- surface, 
in which the multipoles are defined. Another one is a dipole condition, that is concerned 
with the center-of-field and mean-position of the particle. By specifying the dipole in a 
certain manner, we obtained the Papapetrou equations as the equations of motion. Another 
specification led us to more simple representation of the equations of motion for the spinning 
particle. Our equations allow us to understand the motion of the spinning particle in more 
simple way. The deviation from a geodesic can be given by only a coupling term of the spin 
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of the particle and the curvature of background spacetime, and the spin evolves by parallel 
transportation along the world line of the particle. 

The simpler equations bring us the following advantages. In the case of the Papapetrou 
equations, the 4-velocity of the particle is different from its linear momentum in general. 
Therefore, some researchers have imposed the pole-dipole approximation on the equations 
of motion. This approximation is based on the assumption that quadratic terms of the 
spin have the same order of quadrupole moment. However, if we consider a peculiar case 
that an object has spin but no quadrupole moment, for example a spinning rigid body, the 
pole-dipole approximation is no longer valid. Nevertheless, our equations are valid because 
the 4-velocity is parallel to the linear momentum in our equation. 
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APPENDIX A: PROOFS FOR TULCZYJEW'S THEOREMS 



1. Proof for lemma ^ 



Let us focus on the following, 



V p {a^-{r)u p V) 







dxP 



[r)u p V) + T p pa a pu - u p V + ■■■ + Y p pa a pv -u a V 



d , „„ Da pv " 

- ( „r-v ) + —v. 



(Al) 



Note that a M1 "(r) and u p (t) are not function of x and they pass through the partial deriva- 
tives by x. Additionally, 4- and are commutate each other, because r is merely a 
parameter and has no relation to x. Therefore, the first term of the equation (jAljl becomes 



V„---V, 



dr 



{a pu -V) 



dr 



d_ 

d^ 



[V Q ---V^ (a^-V)]dr, 



(A2) 



then we can proof the lemma d 
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2. Proof for theorem [2] 



Our proof is based on the mathematical induction. Firstly, suppose that 
V^j, • • • V M1 n /il '"' tfc , can be reduced to the following form, 

k 

V, fe ■ • • V^™* = E V ^ • • • V^" 1 '"^ 1 , (A3) 

1=0 

with the perfectly symmetric tensor t> Ml ' 

Secondly, let us consider a (k + l)-rank tensor n. w '" M * +1 . We can suppose that it is sym- 
metry with the indices from fix to due to the assumption of the mathematical induction. 
Notice the following equation, 

& + l V / ' 

where we introduce a new notation for anti-symmetrization. The indices with an under- 
line in flA4|) means anti-symmetrization after symmetrization, for instance, n^" -^ 1 = 
i (n&vukhk+i _ n (Mfe+i-Mfc)Mi). Using this relation (|A"4j) . we obtain 

V Mfc+1 • • • V^-^^ = V Mfe+1 • • • V w «C"-^» (A5) 

+FTT ([V ^ +1 ' Vw]Vw • • • V ^ + • • • + ^i. V »1Vi • • ■ V m) n^-^W. 

The terms in the parenthesis of the right hand side of the equation (jA5|) can be written into 

the form of the sum of (R ^ or V ■ • ■ VP"^) x n p "' . We can apply the assumption of the 

induction to these terms. Therefore, we have proofed 

fc+i 

1=0 

Next, let us proof the orthogonality of v^ 1 '"^ 1 and w M . Define a projection operator 
P^ v = S^ u + u^u v . This is orthogonal to w M : P^ v u v = 0. We can write an arbitrary tensor 
X^- in the form X^ 1 '" = ■ • ■ I" 1 '- + ■ • • + u pl (u Ul X Ul -) + Redefine by 

the tensor projected by P^ instead of old v 1 * 1 '" 1 * 1 , then we can assure of the orthogonality 
of new u Atl **'W and 



3. Proof for theorem |3] 



Let us denote the two intersection points of the domain D and the world line of the 
particle by zip) and z(q) (p < q). At the end points zip) and z{q), we can choose the values 
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of K vv .. Vm and its derivatives lower rank than n — 1 to vanish for the sake of arbitrariness of 
K Ul ... Um . Then, substituting ()3.5|) for (j3.6J) and integrating it by parts repeatedly, we obtain 
the following equation, 

[ dr{-l) k [V» k ■ ■■V, 1 ^ 1 .., m ]_ (r) / 1 '-^ 1 -- = 0. (A7) 

fc=0 j p 

Because the values of K Vl ... Vm at z(p) and z(q) are fixed, we can not assume the behavior 
of K ul ... Um along the world line of the particle. However, there is still freedom to fix the 
behaviors of K ul ... Um in the direction perpendicular to the world line. Therefore, we can set 
its spatial derivatives to vanish on the world line. In addition, we change the frame to the 
comoving frame of the particle. Consequently, (|A7|) becomes 

f dr(-ir[d in ■ ■ -d h K„ v ..„ m ] x=z{T f h " 4nn = 0. (A8) 
Jp 

Note that all the derivatives are only spatial ones because of u lll v' J ' 1 "' flnUl "' Um = 0, and also 
note that the indices i\ to i n are fully symmetry because of commutability of the partial 
derivatives. Now, there remains the arbitrariness of the value of \di n ■ ■ ■ di x K Vv .. Vm ^ x _ z ^ T y so 
that, v %v " %nUl '" Vm = is concluded in order for (jA8|) to be valid. This means 



in general coordinates. Repeating the same discussion, it is concluded that all coefficients 
v m-vi- m ust be equal to 0. 
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